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Abstract 

For the L2-orthogonal projection Py onto spaces of linear splines over simplicial 
partitions in polyhedral domains in H*^, d > 1, we show that in contrast to the 
one-dimensional case, where ||-Pv||Loo^ioo — 3 independently of the nature of the 
partition, in higher dimensions the Loo-norm of Py cannot be bounded uniformly 
with respect to the partition. This fact is folklore among specialists in finite 
element methods and approximation theory but seemingly has never been formally 
proved. 



1 Introduction 

Variational methods based on piecewise polynomial approximations are a workhorse in 
numerical methods for PDEs and data analysis. In particular, least-squares methods 
lead to the study of the L2-orthogonal projection operator Py : L2{fl) —>■ V onto a 
given spline space V defined on a domain e R''. A question of considerable interest 
is the uniform boundedness of the Lgo-norm 

of Py with respect to families of spline spaces V. If f2 is a bounded interval in IR^, then 
the question has been intensively studied for the family of spaces of smooth splines of 
fixed degree r over arbitrary partitions, where Shadrin [7j has recently established that 

||Py||L^^L^<C(r)<oo, r>l, (1) 

for any partition. This result was known for a long time for small values of r, e.g., 
Ciesielski [3] proved that for linear splines one can take C(l) = 3, while de Boor [2] 
solved the case r < 4. The estimate plays an important role in numerical analysis 
and for the investigation of orthonormal spline systems such as the Franklin system in 
Lp-based scales of function spaces, 1 < p < oo. 



In higher dimensions, the study of the Loo-norm of Py arose mostly in the context of 
obtaining Lp error estimates in the finite element Galerkin method [5l|6], where sufficient 
conditions on the underlying partition and nodal basis {(pi} of a finite element space V 
are formulated under which the norms ||-Py Hloo^-Lcx) bounded by a certain finite 
constant. Interestingly enough, these results suggest that such conditions on partitions 
resp. finite element type are essential for obtaining uniform bounds but formal proof of 
their necessity was not given. 

Similarly, in the theory of multivariate splines final results on the uniform bound- 
edness of H-Pi/lkcxj^ioo could not be localized. Recently, Ciesielski asked about the 
extension of his result for linear splines [3] to the higher-dimensional case, and the unan- 
imous opinion of the audience was that in higher dimensions a similar result cannot 
hold. However, other than a vague reference to unpublished work by A. A. Privalov, no 
concrete proof could be found. 

It is the aim of this note to provide an elementary example of triangulations 7j of a 
square f2 C IR^ into 0(J) triangles for which 

\\PviTj)\\L^^L^ > J, J>1, (2) 

where V{Tj) is the space of linear C° splines (or finite element functions) on Tj; see 
Theorem [T] below. As J — oo, the triangulations Tj will not satisfy the minimum angle 
condition, which is natural since for these types of triangulations a uniform bound can 
easily be established. However, they satisfy the maximum angle condition, [Ij, and do 
not possess vertices of high valence. The example can easily be extended to d > 3, and 
implies that C(l) = oo for all > 1. 

2 Notation and Result 

We concentrate on d = 2, the case d > 3 is mentioned in Section [31 Let C IR^ be a 
bounded polygonal domain equipped with a finite triangulation T into non-degenerate 
closed triangles A satisfying the usual regularity condition that two different triangles 
may intersect at a common vertex resp. edge only. The set of vertices of T is denoted 
by Vr- Let |£^| denote the Lebesgue measure of a measurable set E C R^. By < a^- < 
ar < TT we denote the minimal and maximal interior angle of all triangles in T. 

Let V{T) denote the linear space of all continuous functions g whose restriction 
to any of the triangles A G T is a linear polynomial. Any g G V{T) has a unique 
representation of the form 

g = Y: 9{P)<Pp. (3) 
PeVr 

where the Courant hat functions </)p G V{T) are characterized by the conditions (t)p{P) = 
1, and (t>p{Q) = 0, where Q ^ P is any of the remaining vertices of T. Thus, dim V{T) = 
#Vr, and 

supp (j)p = fip := UAeT: pgaA. 

The set Qp corresponds to the 1-ring neighborhood of P in T, and we denote by Vp = 
{Q E Qp nV : Q ^ P} the set of all neighboring vertices to P. 
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The L2-orthogonal projection of a function / G L2{fl) onto V{T) is given by the 
unique g := Pv{T)f ^ such that 

{f-gAp) = WPeVr. 

Here and in the sequel, (■, ■) stands for the -^2(^2) inner product. Using with unknown 
nodal values xp = g{P) as ansatz, this is equivalent to the linear system 

$:(0Q,0p)xQ = (/,0p) VPGV. (4) 
Qev 

A simple calculation shows that 

(0p,0p) = 



Similarly, if Q G Vp is a neighbor of P then 

(0Q,0P)= 

Agr:P,QeA -"-^ 

in all other cases we have (0q, 0p) = 0. This shows in particular that 

(0p,0p)= 5](0Q,0p) = (l,0p)/2. 

I.e., if we normalize in by {(f)p,(j)p) then (jl]) turns into a linear system 

Ax = b, X := (XQ : Q G Vr)^, & := {hp = //' '^J^ : P G Vr)^, (5) 
where ||6||oo < 2||/||l^, and A := (apg) satisfies 

f 1, Q = P, 

y 0, otherwise, 

and 

1= E «PQ= E VPG Vt. (7) 

Q^p QeVp 

Thus, the matrix A is only weakly diagonally dominant, and not strictly diagonally 
dominant as in the one- dimensional case. Otherwise, we could estimate ||y4~^||oo in a 
trivial way, and use the inequality 

\\Pv{T)\\l^^l^ < 2p"ioo = 2 niax \\y\\^, (8) 

which follows from the above, in conjunction with the obvious equality ||Pv'(r)/||Loo — 
1 1 x I loo- Let us mention without proof that ([H]) implies the following partial result. 
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Proposition 1 If for any two neighboring vertices P ^ Q from Vr we have apq > cq > 
0, then 

\\Pvir)\\L^^L^<{l + '2co)c^\ 

For triangulations satisfying ttie minimum angle condition uniformly, i.e., a-^- > > 0; 
this result is applicable with a cq determined solely by cto, and thus covers the bounds 
considered in the finite element literature [5], [6]. The main result of this note is the 
following 

Theorem 1 For any J > I, there is a triangulation Tj of a square into 8J + 4 triangles 
such that the norm of the L2-projector Pv{rj) satisfies 

\\Pv{rj)\\L^^L^ > 2 J. 
Thus, for spatial dimension d = 2 we have 

C(l) := sup 1 1 Py(r) 1 1 Loo ^Loo = oo- 

T 

We conjecture that in terms of the number of triangles our result is asymptotically 
sharp for bounded polygonal domains in IR^, i.e., 

sup \\Pv{T)\\l^^l^-N, iV^oo. (9) 

T: #T<N 

That is a square is not crucial. The examples given below can easily be modified to 
show C(l) = oo for simplicial partitions in higher dimensions as well. 



3 Proof of Theorem [T] 

We will use the following notation. Let Sa = [—a, a]^ be the square of side-length 2a, 
a > 0, with center at the origin. Set fl := Si, with vertices denoted (in clockwise 
direction) by Po.i, i = 1, • • • ,4, and fix the parameter t G (0, 1). The triangulation 7j, 
J > 1, of is obtained by inserting the squares St, St2,. . . , StJ, whose vertices will be 
denoted similarly by Pj^, i = 1, . . . , 4, j = 1, . . . , J, placing an additional vertex Pj+i 
at the origin, connecting Pj+i by straight lines with the 4 vertices Po,j of Si, and finally 
subdividing each of the remaining trapezoidal regions Pj_i^j_iPj- j_iPj^jPj_i^j into two 
triangles by connecting Pj_i i_i with Pji in a consistent way. The outer rings of the 
resulting triangulation are shown in Fig. [T]a). 

Let the function / G L^oi^) with ||/||loo = 1 be defined as follows: 



fix) = i-iy, xeQr- 



StJ-l\St3, j — 1, . . . , J, 

StJ, j = j+i. 



We use the same notation as above, and consider the linear system ([5]) corresponding to 
the L2-orthogonal projection g = Pv{rj)f of this / onto V(Tj). Because of uniqueness of 
orthogonal projections and the rotational symmetry of Tj and /, the entries of the vector 
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Figure 1: Triangulation Tj for t = 0.3 (left), and typical fip^. . (right) 

X that represent the nodal values of g corresponding to the vertices Pj^i, i = 1, ... ,4, 
are equal, and will be denoted hj Xj, j = 0, . . . , J (the value at the origin is denoted 
by xj+i). Moreover, the 4 equations in ([5]) corresponding to the 4 vertices of a square 
St^, j = 0, . . . , J, can be replaced by one, thus turning the original system Ax = b of 
dimension 4 J + 5 into a reduced tridiagonal system Ax = b of dimension J + 2 for the 
vector X = (xq, . . . , x j+i)"^. Since the equations in the system ([5]) are invariant under 
affine transformations, and because of the definition of Tj via squares shrinking at a 
fixed geometric rate, it is easy to see that, with the exception of the first and last two, 
all equations in Ax = b have the same form: 

axj^i + /3xj + 'yxj+i = {-lyS, j = 1, . . . , J - 1. (10) 

The coefficients can be found from the triangular neighborhood fip - of any of the Pj^i, 
j = 1, . . . , J — 1 (see Fig. [T]for an illustration), and the definitions leading to (ISl). We 
do not need their exact coefficient expressions, just their limit behavior as t 0, i.e., 
we will be looking for the entries of A := limj^o^ and b := limt^ob. Indeed, since for 
t the whole is essentially covered by the single triangle with vertices Pj^i, Pj-i^i-i, 
Pj-i,i, and since f{x) = (—1)-' on the latter, we have 

a = l + 0(t), (3=l + 0{t), 7 = 0(^2), (5 = 2 + 0(t). 

From (ITO!) . we obtain in the limit t — the equations 

Xj-i + Xj = 2{-iy, j = i,...,j-i, 

where Xj = lim^^oo Xj (the existence of these limits follows from the invertibility of the 
limit matrix A, see below). Similar considerations for the first and the last two equations 
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of Ax = h yield the remaining three equations of Ax = b as follows: 

3. 1 



Xo + -xi = -2, Xj_i + Xj = 2{-iy , j = J,J + l. 



The resulting matrix A is obviously invertible. After finding Xq — Xi — — 1 from the 
first two equations, forward substitution gives Xj = (2j — 1)(— I)-', j = 2, . . . , J+ 1. This 
implies that for any e > one can find a sufficiently small t > such that 

\\x\\oo > \\x\\oo - e = 2J + 1 - e. 

This proves Theorem [H Note that the above reasoning does not work for type-I trian- 
gulations of a square obtained from a non-uniform rectangular tensor-product partition. 

We conclude with the straightforward extension of the above example to arbitrary 
d > 2. Let e™ denote the m-th unit coordinate vector in R'^, m = 1, . . . , rf. As we 
take the convex polyhedral domain with vertices 

Po,l = e^ + e^ Po,2 = e^-e^ Po,3 = -e' - e^ Po,4 = -e' + e^ 

and = e™, m = 3, . . . , (i. For d = 3, this domain is a pyramid with square base in 
the xy-plane, and tip on the z-axis. 

A suitable simplicial partition of O is obtained as follows. The base square with 
vertices Po,i; i = l;---;4, is triangulated into Tj which depends on the parameters 
< t < 1 and J as described for d = 2. The resulting triangulation (now embedded into 
R'^) and its vertices are again denoted by 7} resp. by Pj j, i = 1, . . . , 4, j = 0, . . . , J, 
and Pj+i. Then each simplex in the associated simplicial partition Vj of Q is generated 
by the d — 2 vertices P^, m = 3, . . . ,d, and the three vertices of a triangle in Tj. The 
latter is called base triangle of the associated simplex. Obviously, the ci-dimensional 
volume of each simplex is proportional to the 2-dimensional area of its base triangle, 
with proportionality constant 2/d\. To obtain a suitable function / G Loo(f^) with 
ll/IUoo = 1 W6 prescribe values ±1 on the simplices by inheritance from the values on 
the base triangles of the above 2-dimensional /. 

From the symmetry properties of / and Vj, it is obvious that the L2-orthogonal 
projection Pv(Vj)f of / onto the linear spline space V(Vj) is characterized by its value 
xj+i at the origin Pj+i, by values Xj taken at the vertices Pj^i, i = 1,...,4, where 
j = 0, . . . , J, and a common value x' taken at the remaining vertices P^, m = 3, . . . ,d. 
To estimate these values which, in complete analogy to the 2-dimensional case, are 
represented by the solution vector x of a certain linear system Ax = b (now of dimension 
J + 3), we need the limit version Ax = 6 of this linear system for t —>■ 0. We spare the 
reader the elementary calculations, and state it without proof: 

d-2 , , ,.id + 2 
Xj+Xj-i^ = (-1)^—, J = 1, •••,</+ 1, 

3. 1. rf-2., d + 2 

-Xq H Xi H X = , 

2 2 2 2 ' 

1. d-3,^, d + 2 
Xq + -Xi + (1 H — )x = — . 
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Prom this system one easily concludes that ||£||oo > cJ, and consequently ||-Py(-Pj)/||Loo — 
cJ for a small enough t > which implies the desired result. The lower bound cJ 
obtained does not seem to accurately reflect the possible growth of the projector norms 
in Loo as a function of the number of simplices, for d > 3 we would rather expect an 
exponential rate. 
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